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Abstract: The reactivities of the ada- 
mantane-like heteronuclear vanadium- 
phosphorus oxygen cluster ions 
[V JC P 4 _ JC O 10 ] ,+ (x = 0, 2-4) towards hy- 
drocarbons strongly depend on the V/P 
ratio of the clusters. Possible mecha- 
nisms for the gas-phase reactions of 
these heteronuclear cations with 
ethene and ethane have been elucidat- 
ed by means of DFT-based calcula- 



tions; homolytic C— H bond activation 
constitutes the initial step, and for all 
systems the P— O' unit of the clusters 
serves as the reactive site. More corn- 
Keywords: C— H activation • cluster 
compounds • density functional cal- 
culations • gas-phase reactions • re- 
action mechanisms 



plex oxidation processes, such as 
oxygen-atom transfer to, or oxidative 
dehydrogenation of the hydrocarbons 
require the presence of a vanadium 
atom to provide the electronic prereq- 
uisites which are necessary to bring 
about the 2e" reduction of the cationic 
clusters. 



Introduction 

Selective oxidation of hydrocarbons continues to constitute 
one of the major challenges in contemporary chemistry to 
solve global problems, such as an environmentally benign 
and economically feasible conversion of natural gas into 
value-added products. [1] In this context, oxygen-based cata- 
lysts, especially metal oxides, play a particular role as power- 
ful oxidation reagents. [2] Despite the apparently large 
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number of different pathways and mechanistic scenarios, 
formation of the various oxidation products of hydrocarbons 
by metal oxides can be classified in terms of three general 
reaction types (Scheme 1): i) hydrogen-atom transfer (HAT) 
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Scheme 1. Schematic description of three basic oxidation processes of hy- 
drocarbons by metal oxides [MO] , all commencing with C— H bond acti- 
vation. 



from the hydrocarbon to the metal oxide to bring about oxi- 
dative coupling, ii) oxygen-atom transfer (OAT) to the or- 
ganic substrate, and iii) oxidative dehydrogenation (ODH) 
of the latter. While there exist numerous effective homoge- 
neous and heterogeneous catalysts which cover this rather 
broad spectrum of oxidation reactivity, it is no exaggeration 
to state the paucity of detailed knowledge about the intrin- 
sic properties of many of the catalysts which, after all, con- 
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trol also the chemoselectivity of the three competing oxida- 
tion processes depicted in Scheme l. [3] 

To further improve the understanding on how such com- 
plex systems, for example, surfaces in heterogeneous cataly- 
sis may operate, often model systems have been studied as 
they can provide insights into what happens at the molecu- 
lar level. [4] One approach in the elucidation of mechanistic 
aspects employs reactivity studies which are conducted 
under near single-collision conditions in the gas phase. In 
fact, in the last two decades numerous state-of-the-art, 
mass-spectrometry based experiments with small metal-oxo 
clusters in combination with computational studies have re- 
vealed important details of the elementary steps for many 
chemical reactions. [5] For example, hydrogen-atom abstrac- 
tion from CH 4 to generate CH 3 ' is viewed as the decisive 
step in the oxidative dehydrogenation and dimerization of 
methane, [6] and many studies suggest the crucial role of 
oxygen-centered radicals to bring about homolytic C— H 
bond scission by (doped) metal oxides. [2b 7] The perhaps 
most compelling experimental evidence for this concept has 
been provided by the reactions of free aluminum-oxide clus- 
ter ions with methane : [8] aluminum-oxide cluster cations 
with an odd number of aluminum atoms, thus lacking an 
oxygen-centered radical, do not react with methane at ambi- 
ent conditions, in spite of a substantial thermodynamic driv- 
ing force; in distinct contrast, oligomeric [(A1 2 0 3 )J' + (x = 3, 
4, 5) cluster ions possessing an even number of aluminum 
atoms and an oxygen-centered radical are remarkably reac- 
tive towards methane even at room temperature. 

In addition to the large number of homonuclear metal 
and non-metal oxides investigated^ 8 ' 9] more recently hetero- 
nuclear oxide clusters have been brought into focus, for ex- 
ample, [M x V y O z ] +l - (x+y = 2, 3, 4; z = 3-10), [10] [Ce x V y O z ] + 
(x+y = 2, 3; z = 4, 5, 6),^ or [(V 2 0 5 ) x (Si0 2 )J +/ - (x = l, 2; 
y = l-4), [12] respectively. The considerable interest in study- 
ing these systems as model catalysts can be derived from 
their chemically close resemblance with practical heteroge- 
neous catalysts. Bulk metal-oxide catalysts are often com- 
posed of multiple components, that is, mixed-metal oxo- 
frameworks, or employ metal oxides as non-innocent sup- 
port materials. [13] The transformation of ^-butane to maleic 
anhydride by the so-called VPO-catalysts (Scheme 2), may 

o 

VPO-catalyst_ rrA 



Scheme 2. Industrial synthesis of maleic anhydride from ^-butane. 

serve as a well-known example for such a multi-component 
system in heterogeneous catalysis; this reaction involves the 
abstraction of eight hydrogen atoms from as well as the 
transfer of three oxygen atoms to C 4 H 10 — complexity at its 
best. [14 ^ Since the first commercial use of VPO catalysts for 
the synthesis of maleic anhydride in the early 80s, many ex- 
perimental and theoretical studies have been carried out to 
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investigate this remarkable transformation. [15] However, the 
underlying reaction mechanisms are still only poorly under- 
stood and many conflicting hypotheses exist. For example, it 
was conjectured that structurally ill-defined phosphate spe- 
cies work as catalytically innocent, chemically inert linkers 
between the active vanadium-oxide sites in these mixed 
metal-oxide phosphates ; [14a] however, this assumption was 
questioned by a gas-phase study on the open-shell oxide 
cation [P 4 O 10 ]' + , in which the efficient activation even of 
methane, the most inert of all hydrocarbons, at ambient con- 
ditions by this metal-free oxide was demonstrated. [9d] Obvi- 
ously, the central question about the actual nature of the 
active site in heteronuclear oxo-clusters and the particular 
role of the cluster components was once more raised. 

Additional gas-phase studies supported the idea that also 
phosphorus sites may be active by demonstrating an in- 
creased reactivity of a polynuclear metal-free cluster to- 
wards various small hydrocarbons, for example, ethane, pro- 
pane, and ethene. [16] This phosphorus-oxide system is elec- 
tronically and structurally closely related to the isostructural 
transition-metal oxide cluster [V 4 O 10 ]' + , which had been 
thoroughly investigated earlier; it seems to be an appropri- 
ate model for surface-mediated C— H bond activation proc- 
esses.^ Common to both cluster-oxide cations is the exis- 
tence of a terminal oxygen-centered radical which proved to 
be crucial for the reactivity observed towards methane and 
higher hydrocarbons. [9ad l618] While [P 4 O 10 ] ,+ and [V 4 O 10 ] ,+ 
bring about thermal hydrogen-atom abstraction from meth- 
ane, the reactivity patterns with ethane, propane and ethene 
differ quite remarkably. For example, [P 4 O 10 ]' + reacts with 
all substrates predominantly by homolytic C— H bond cleav- 
age (HAT), while [V 4 O 10 ]' + gives rise to OAT to the higher 
hydrocarbons as well as to ODH. [19 ^ These differences in ox- 
idation reactivities suggested a more systematic investigation 
at a molecular level, that is, the studies of mixed vanadium- 
phosphorus oxygen cluster ions [V x P 4 _ x O 10 ]* + (x = l-3). [20] 
While only the mixed cluster ions [V 2 P 2 O 10 ]' + and 
[V 3 PO 10 ]' + could be generated in the experiment, density 
functional calculations revealed for all three mixed clusters 
[VP 3 O 10 ] ,+ , [V 2 P 2 O 10 ] ,+ , and [V 3 PO 10 ] ,+ a phosphorus-bound 
oxygen-centered radical as the most stable isomer on the 
potential-energy surfaces (PESs). Thus, all C— H bond acti- 
vation processes are initiated solely by the terminal phos- 
phorus-bound oxygen atom. Not entirely unexpected with 
regard to the activation of methane by both oxides [P 4 O 10 ]' + 
and [V 4 O 10 ]' + , the reactivities of [V 3 PO 10 ]* + and [V 2 P 2 O 10 ]' + 
towards methane are similar in that hydrogen-atom abstrac- 
tion is observed with comparable efficiencies. [9ad ' 20b] Howev- 
er, the reactions of [V 3 PO 10 ]' + and [V 2 P 2 O 10 ]' + with 
C 2 hydrocarbons give rise to new product distributions for 
the couples [V 3 PO 10 ]' + /C 2 H x and [V^O^^/QH, (jc = 4, 6). 
While homolytic C— H bond scission, according to Equation- 
s (la) and (2a), and oxidative dehydrogenation, [Eqs. (lc) 
and (2c)] , are observed, oxygen-atom transfer does not take 
place [Eqs. (lb) and (2b)] . 

The questions raised by the experimental data (Table 1) 
can be summarized as follows : i) What are the reasonable 



3018 



www.chemeurj .org 



© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim 



Chem. Eur. J. 2013, 19, 3017-3028 



Gas-Phase Reactions of Vanadium-Phosphorus Oxide Clusters with C 2 H X 



FULL PAPER 



[V 3 PO 10 r + C 2 H X 









— # — ► 



[V 3 P0 9 (OH)] + + C 2 H x .i 



(la) 



-H ► [V 3 P0 9 ]* + + C 2 H x O (lb) 



[V 2 P 2 O 10 ]' + + C 2 H X 





— # — ► 







[V 3 PO 10 H 2 p + C 2 H X _ 2 

[V 2 P 2 0 9 (OH)] + + C 2 H x .r 



-n — ► [v 2 p 2 o 9 ]* + + c 2 h x o 



(lc) 
(2a) 

(2b) 



[V 2 P 2 O 10 H 2 ]" + + C 2 H X _ 2 (2c) 



Table 1. Branching ratio of the product distributions in the reactions of 
[V^P 4 x O 10 ]' + (x = 0, 2-4) (OAT = oxygen-atom transfer, HAT = hydro- 
gen-atom transfer, ODH = oxidative dehydrogenation). [1618b ' 20b] 
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[a] According to our experimental results. [19] [b] Minor contributions of 
single-electron transfer to produce [C 2 H 4 ]' + have not been taken into ac- 
count. 



As shown in Table 2, BDE(V— O t ) increases when phos- 
phorus atoms are present in the cluster: In [V 4 O 10 ]' + (ly 4 )» 
BDE(V— O f ) amounts to only 280kJmol _1 , which is 
76kJmol" 1 and 71 kJmol -1 lower compared to the values 
obtained for the [V 3 PO 10 ] ,+ (l VjP ) and [V 2 P 2 O 10 ] ,+ (1 V2 p 2 ) 
clusters, respectively.^ In contrast, the BDEs of the termi- 
nal P— O bonds for all phosphorus-containing clusters inves- 
tigated are not affected by the vanadium/phosphorus ratio 
and lie around 400 kJmol" 1 . These features can be explained 
by the fact that the most stable structures of all phosphorus- 
containing clusters possess a P— O' bond; a V— O' bond is 
only present in [V 4 O 10 ]' + . Thus, energy in addition to BDE- 
(V— O f ) = 280 kJmol" 1 as determined for [V 4 O 10 ]' + is neces- 
sary for [V 3 PO 10 ]* + and [V 2 P 2 O 10 ]* + to cleave the inert V=0 
bond. [24] While for the heteronuclear clusters it is energeti- 
cally favored to remove an electron from the P=0 bond, 
generating a P— O' instead of a V— O' unit, the binding 
energy of the P— O' versus the V— O' bond is higher. An anal- 
ysis of the nature of the singly occupied molecular orbitals 
(SOMOs) of the dissociation products [V x P 4 _ x 0 9 ] ,+ (x=2, 3) 
is quite helpful to understand the origin of these thermo- 
chemical differences; these SOMOs are utilized to form a 
single bond to oxygen in [V x P 4 _ JC O 10 ]' + . As exemplified in 
Figure 1 for [V 2 P 2 O 10 ]' + , loss of a terminal O atom from va- 
nadium leads to a vanadium-centered radical having a 



reaction mechanisms for the various VPO-cluster 
ions in their thermal reactions with ethene and 
ethane? ii) Why does only the [V 4 O 10 ]' + /C 2 H 4 - 
couple undergo exclusively oxygen-atom-transfer, 
but as soon as one phosphorus atom is present in 
the tetranuclear oxide cluster, oxygen-atom trans- 
fer to the hydrocarbon is no longer observed? 
iii) Why does [P 4 O 10 ]' + not permit the oxidative de- 
ny drogenation of the C 2 -hydrocarbons ? 

Herein, we aim at providing a DFT-based explanation for 
the different reactivities exhibited by [V x P 4 _ x O 10 ]* + (x = 0, 2- 
4), or to put it bluntly "why composition counts! 



Table 2. BDEs (kJmol 1 ) of the terminal V-O and P-O bonds in [V X P 4 x O w ]' + (x = 0, 
2-4), calculated at the B3LYP/aug-cc-pVTZ//B3LYP/TZVP level of theory. 
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Results and Discussion 

Intrinsic properties of the [\ x Y 4 _ x O 10 \' + clusters (x=0> 2-4): 

Before discussing the relative energies of the intermediates 
and transition structures as well as the mechanistic details of 
the three reaction channels for the individual [V JC P 4 _ JC O 10 ]' + 
/C 2 H 4 and [V JC P 4 _ J A 0 ]' + /C 2 H 6 systems (* = 0, 2-4), some of 
the intrinsic properties of the individual clusters are ad- 
dressed; this may be helpful in the discussion of the experi- 
mentally observed reactivity patterns. [21] For example, to ex- 
plain the preference for OAT observed in the reaction of 
[V 4 O 10 ]' + with ethene, and its complete absence for the 
phosphorus-containing systems, a closer look at the bond- 
dissociation energies (BDEs) for the terminal V— O f and 
P— O f bonds, respectively, is instructive (O t = terminal 
oxygen). [22] 



P=0 



v=o 



'4x 



Figure 1. Two possible DFT-derived [V 2 P 2 0 9 ]' + isomers resulting from 
dissociation of a P— O' or V— O* bond of [V 2 P 2 O 10 ]' + (green V, yellow P, 
red O). The blue isosurfaces indicate the spin density within the respec- 
tive cluster. 



d z i orbital as a SOMO. In contrast, cleavage of the P— O/ 
bond in [V 2 P2O 10 ]' + gives rise to a phosphorus-centered 
SOMO possessing a s/p character (Figure 1); the latter is in 
comparison to the d orbital of the vanadium atom, more ap- 
propriate to form a single bond to the oxygen atom due to 
the better overlap with its sp* orbital ; similar properties of 
the SOMO were obtained for all P-containing cluster ions. 
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The thermochemistry of the competing HAT and 
ODH reaction channels depends on the relative en- 
ergies of the [V JC P 4 _ JC O 10 H] + and [V JC P 4 _ x O 10 H 2 ]' + 
product ions, respectively. While the first hydrogen- 
atom abstraction results in the formation of closed- 
shell product ions concomitant with spin transfer to 
the emerging hydrocarbon radicals, respectively, 
open-shell systems are regenerated in the course of 
the second hydrogen-atom transfer under reduction 
of a vanadium or phosphorus atom, respectively, to 
the formal oxidation state +IV. As will be discussed 
later, the reduction to the oxidation state +IV is 
disfavored for phosphorus. Thus, [P 4 O 10 ]' + gives rise 
to the closed-shell HAT product [P 4 0 9 (OH)] + ; for- 
mation of the open-shell product [P 4 O 10 H 2 ]' + via 
ODH is not possible due to thermochemical con- 
straints. In contrast, the presence of redox-active 
vanadium opens up the latter reaction channel. The 
relative stabilities of the open- versus the closed- 
shell cations are indicated by the adiabatic ioniza- 
tion energies (IEs) of the respective neutral clusters 
(Table 3). The formation of open-shell cluster cations, that 
is, [V x P 4 _ x O 10 ]' + (1), [V x P 4 _ x 0 9 r, and [V x P 4 _ x O 10 H 2 ]' + , is 
rather energy demanding in the case of the homonuclear 
phosphorus systems (x = 0, Table 3) ; in contrast, the ioniza- 
tion energy required to produce the closed-shell cluster 
cation [P 4 O 10 H] + is significantly lower (7.2 eV) than the IEs 
of the respective vanadium-containing systems (>8.5 eV). 



Table 3. Adiabatic ionization energies (in eV) of V x Y 4 _ x O 10 , V JC P 4 _ x 0 9 , 
V X P 4 x O 10 H and V X P 4 x O 10 H 2 as calculated at the B3LYP/aug-cc-pVTZ// 
B3LYP/TZVP level of theory. 



Table 4. Relative energies and free energies (in parenthesis), given in kJmol -1 , of the 
transition states and intermediates for oxygen-atom transfer (OAT), hydrogen-atom 
transfer (HAT), and oxidative dehydrogenation (ODH) for the reactions of 



[V X P 4 x O l0 Y (x = 0, 2-4) with C 2 H 4 , 
pVTZ//B3LYP/TZVP level of theory. 



respectively, calculated at the B3LYP/aug-cc- 
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[a] The ground state of the neutral cluster is a triplet, while the ground 
states of all other neutral clusters in the respective row are closed-shell 
singlets. 



Ion/molecule reactions with C 2 H 4 : Having discussed some of 
the thermochemical properties of the educt and product 
ions, next we will outline the mechanistic details for the re- 
actions of the [V 2 P 2 O 10 ]' + /ethene system; this may serve as a 
representative example for all the VPO-cluster ions investi- 
gated (the Cartesian coordinates of all structures can be 
found in the Supporting Information, electronic energies 
and the relative free energies are summarized in Table 4), 
and the corresponding potential-energy profile is shown in 
Figure 2. 

In the first step, a carbon atom of ethene coordinates in a 
radicophilic fashion to the terminal phosphorus-bound 
oxygen atom of the cluster to form intermediate 2 in which 
the spin density is mainly located at the terminal carbon 
atom for all of the vanadium-containing intermediates 2 V4 , 



2 V3 F> and 2 V2 p 2 (0.866, 0.893, and 0.866, respectively); in con- 
trast, the spin is distributed over the two carbon atoms in 2p 4 
(0.310 at the proximal and 0.608 at the terminal carbon 
atom). A comparison of the C— O, C— C, and P— O bond 
lengths in 2 for the individual clusters also reveals a notable 
effect due to the presence of vanadium. For all systems in- 
vestigated, r(C— O) in 2 is longer compared to a C— O bond 
in for example, ethanol (1.431 A); however, while a smooth 
increase of r(C— O) is observed in going from 2 V4 via 2y 3 p to 
2 V2 p 2 (1.517 A, 1.583 A, and 1.619 A, respectively), r(C-O) 
in the pure phosphorus intermediate 2p 4 is with 2.021 A 
much more elongated. Though less pronounced, decreasing 
trends in the same direction are observed for the C— C and 
P— O bond lengths of the four systems. The addition of 
ethene to the cluster shortens the P— O bond from 1.599 A 
in ly 3 p to 1.539 A in 2 V3 p and from 1.593 A in 1 V2 p 2 to 
1.530 A in 2 V2 p 2 , respectively; the C— C bonds of 2 are length- 
ened from 1.326 A in free C 2 H 4 to 1.461 A in 2 V4 , 1.455 A in 
2y 3 p and 1.449 A in 2 V2 p 2 , respectively. In contrast, the short- 
ening of the P— O bond is the highest for the [P 4 O 10 ] + /C 2 H 4 
cluster (from 1.569 A in J_p4 to 1.493 A in 2p 4 ) while the elon- 
gation of the C— C bond is the smallest for this particular 
system (from 1.326 A in free ethene to 1.408 A, respective- 
ly); thus, the P— O and the C— C bonds are the shortest for 
the vanadium-free intermediate 2p 4 . 

Starting from intermediate 2, two different reactions 
channels A (OAT) and B (HAT and ODH) are accessible 
(Figure 2). In pathway A, which has already been identified 
by Castleman and co-workers for the reaction of [V 4 O 10 ]' + 
with C 2 H 4 , [18b] a hydrogen atom from the oxygen-bound 
methylene unit undergoes a 1,2 migration to the terminal 
carbon atom via transition structure TS2-3 to form inter- 
mediate 3. The binding situation of the newly generated ace- 
taldehyde building block to the cationic cluster fragment in 
3 differs significantly for the pure phosphorus cluster as 
compared with the vanadium-containing systems: While the 
C-O bond lengths in 3 V4 , 3 V3 p and 3 V2 p 2 (1.242 A, 1.262 A, 
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Figure 2. Potential-energy surface (PES) for the reactions of [V 2 P 2 O 10 ]' + with C 2 H 4 , calculated at the B3LYP/aug-cc-pVTZ//B3LYP/TZVP level of theory 
(green V, yellow P, red O, gray C, white H). The electronic energies and relative Gibbs free energies (in parenthesis) are given in kJmol 1 and corrected 
for unsealed zero-point energy contributions. 



and 1.266 A, respectively) are similar to that of r(C— O) in 
free acetaldehyde (1.205 A), r(C-O) in 3p 4 is with 1.472 A 
relatively long; these structural characteristics are inversely 
mirrored in the shortening of the P— O bond lengths to 
1.745 A, 1.733 A, and 1.526 A in 3y 3 p, 3y 2 p 2 , and 3p 4 , respec- 
tively. Remarkable differences among the clusters investigat- 
ed also concern the spin distributions of the intermediates 
and transition structures for the intramolecular hydrogen 
migration, that is, 2^TS2-3^3. In line with avoiding a high 
spin density at a P atom and the energetic disadvantage of 
an associated reduction to the formal oxidation state P +v ^ 
P +IV as described above, the spin in 3p 4 is located at the 
C2 unit, that is, at the proximal carbon atom of the ligand 
(Scheme 3) ; here, the transfer of the unpaired electron takes 
place only during the reductive elimination of acetaldehyde, 
3p 4 ^4p 4 . The bonding situation is different for the vanadi- 
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Scheme 3. Schematic description of the process 2^TS2-3^3 for the ho- 
monuclear [X 4 O 10 ]' + /C 2 H4 couples (X = P, V). 



um-containing couples [V JC P 4 _ x 0 1 o]' + /C 2 H4 (x = 2-4): In 
TS2-3 the spin density has already been shifted to a vanadi- 
um atom of the cluster; thus, the actual hydrogen transfer 
2^TS2-3— >3 corresponds rather to a 1,2 hydride shift within 
a cationic, closed-shell C2 unit (the C2 units in transition 
structures TS2-3 possess a positive charge for all systems in- 
vestigated), as exemplified for TS2-3 V4 in Scheme 3; this 
process is energetically favored as compared to the energy- 
demanding 1,2 H-atom migration within the radical-cation 
C2 unit in TS2-3p 4 . 

The alternative to the OAT pathway shown in Figure 2, 
that is, the transfer of an oxygen atom from a vanadium- 
bound oxygen of [V 2 P2O 10 ]' + is thermodynamically favored; 
the same holds true for [V 3 PO 10 ]' + (see Table 2). Formation 
of the resulting products 4 V2 p 2 — VO and 4 V3 p— VO are more 
exothermic (exergonic) by —45 (—43) kJmol" 1 and —41 
(—41) kJmol" 1 , respectively, than the corresponding product 
ions of the OAT mechanism depicted in Figure 2. However, 
a rather energy demanding intramolecular spin-density 
transfer from the P=0 to the V=0 site is required, [20a] and 
the associated adduct complexes 2 V2 p 2 — VO and 2 V3 p — VO, in 
which ethene is coordinated to the VO side of the cluster, 
are 59 (56) kJmol" 1 (2 V2 p-VO) and 59 (54) kJmol" 1 
(2 V3 p-VO) less stable than the P=0-bound isomers 2 V2 p 2 and 
2y 3 p respectively. The transition structures TS2-3 V2 p 2 — VO 
and TS2-3 V3 p-VO are with -141 (-101) kJmol 1 and -120 
(—81) kJmol" 1 (relative to the reactants, respectively) com- 
parable in energy to the barriers of the OAT pathway A de- 
picted in Figure 2. Thus, both OAT variants are accessible at 
thermal conditions for the mixed clusters; note, however, 
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that a kinetic barrier caused by the intramolecular spin-den- 
sity transfer might be involved in the initial coordination of 
C 2 H 4 to a non-radical V=0 site of the cluster, as reported 
for the reaction of [V 3 PO 10 ]' + with methane. [20a] Regarding 
the electronic structures of transition states TS2-3 V2 p 2 — VO 
and TS2-3 V3 p-VO, the depicted situation for TS2-3 V4 in 
Scheme 3 holds true also for these systems ; the process cor- 
responds rather to a 1,2 hydride shift within the cationic hy- 
drocarbon fragment as the spin is mainly localized at the 
proximal vanadium atom. 

Pathway B results either in HAT or ODH and commences 
with a hydrogen-atom transfer from the oxygen-bound 
methylene unit in 2 to this very oxygen atom (TS2-5) thus 
generating a oxo-hydroxo cluster with an ethenyl radical 
ligand; the latter is loosely bound to the cluster via the hy- 
drogen atom of the newly generated OH group. In the 
energy profile, Figure 2, the energy of TS2-5 is slightly lower 
than the energy of the resulting intermediate 5; this is due 
to zero-point energy and free-energy corrections. The fact 
that transition state TS2-5 and intermediate 5 are very close 
in energy is reasonable, since the only structural reorganiza- 
tion in this step corresponds to the widening of the binding 
angle of the C 2 H 3 fragment to the cluster. This HAT-reac- 
tion path has been located for all vanadium-containing clus- 
ters, while TS2-5p 4 could not be located for the [P 4 O 10 ]' + 
/C 2 H 4 -couple on the potential-energy surface (PES). Fur- 
thermore, the nature of intermediate 5p 4 differs remarkably 
from the structures of 5 V2 p 2 , 5 V3 p and 5 V4 . In the latter inter- 
mediates, the C— H bond is already strongly elongated 
(r(C-H) 1.675 A for 5 V2 p 2 , 1.756 A for 5 V3 p and 1.643 A for 
5y 4 ) and the newly formed CKH bond is rather short (r(0- 
H) 1.053 A for 5 V2 p 2 , 1.030 A for 5 V3 p and 1.051 A for 5 V4 ). In 
5p 4 , however, the C— H bond is only slightly elongated 
(r(C-H) = 1.102 A) and the CKH distance amounts to 
1.858 A. Thus, 5p 4 rather corresponds to a hydrogen-bridg- 
ed [25] adduct complex [P 4 O 10 -H-C 2 H 3 ]' + directly generated 
from the educts instead of being formed via 2p 4 ^TS2-5p 4 ^ 
5p 4 ; similar adduct complexes could not be located for the 
vanadium-containing clusters with ethene. According to a 
relaxed scan of the C— H bond, the ethenyl radical C 2 H 3 * is 
liberated from 5p 4 without the involvement of any further in- 
termediates or transition structures. In addition to the struc- 
tural differences of 5p 4 as compared with 5 V2 p 2 , 5 V3 p and 5 V4 , 
also the formation of 5p 4 is with — 224kJmol _1 significantly 
more exothermic than that of the other HAT-intermediates 
(-142kJmol" 1 for 5 V2 p 2 , -129kJmol" 1 for 5 V3 p and 
— 123kJmol _1 for 5 V4 ). In line with these findings, the 
energy necessary for the liberation of the C 2 H 3 ' radical from 
5 to form the hydrogen-atom transfer (HAT) product 6 
amounts to 124 kJmol" 1 for [P 4 0 9 (OH)] + (6p 4 ), compared to 
only 55kJmol 1 for both [V 2 P 2 0 9 (OH)] + (6 V2 p 2 ) and 
[V 4 0 9 (OH)] + (6 V4 ), and 47kJmol" 1 for [V 3 P0 9 (OH)] + 
(6 V3 p), respectively. 

Structure 5 V2 p 2 in Figure 2 constitutes to a common inter- 
mediate for both the HAT and the ODH pathways. In com- 
petition with the formation of the HAT product 6, the incip- 
ient C 2 H 3 ' radical can rebind to the oxygen atom of the 



newly formed hydroxyl group thus producing structure 7 V2 p 2 ; 
due to the rather shallow nature of the PES in this region, a 
transition state TS5-7 could not be located on the PES. 
However, the unsuccessful attempts to locate this transition 
state suggest that the barrier for this rebound step is lower 
than the liberation of the hydrocarbon fragment resulting in 
the HAT product 6. Regarding the structure of clusters 7, all 
vanadium-containing cations 7 V2 p 2 , 7 V3 p and 7 V4 possess the 
bonding pattern depicted in Figure 2 with the spin density 
exclusively located at one of the vanadium atoms within the 
cluster skeleton. For the pure phosphorus system instead, 
the rebound step is associated with an opening of the cluster 
cage thus leading to a distorted structure of 7p 4 (see the Sup- 
porting Information). For the vanadium-containing clusters 
the reaction continuous via transfer of the C 2 H 3 fragment to 
a bridging, vanadium-bound oxygen atom. In the case of 
[V 2 P 2 O 10 ]' + , the C 2 H 3 fragment can alternatively be transfer- 
red to the bridging, phosphorus-bound oxygen atom; this 
pathway is described further below. Note also that starting 
from 5 the direct migration of the C 2 H 3 * radical to a bridging 
oxygen atom also results in the formation of 8 and repre- 
sents an alternative mechanistic scenario; however, locating 
the corresponding transition structure was not successful. 
From 8, transfer of a hydrogen atom from the terminal 
CH 2 unit to the adjacent V=0 site, thus generating a V-OH 
group, completes the ODH process. The barrier height of 
the step 8^TS8-9^9 depends strongly on the system, rang- 
ing from 84 and SSklmol" 1 for [V 2 P 2 O 10 ] ,+ and [V 3 PO 10 ] ,+ , 
respectively, to 151 kJmol" 1 for [V 4 O 10 ]' + . For the latter, the 
high intrinsic barrier is due to the relative low energy of in- 
termediate 8y 4 , in which the bond length of bridging u.- 
OCHCH 2 ligand to the V(OH) and VO, units, respectively, 
are quite similar (r((HO) V-OCHCH 2 ) = 2.003 A, 
r(O y V-OCHCH 2 ) = 1.918 A) ; in contrast, the u,-OCHCH 2 
ligand in 8 V2 p 2 and 8 V3 p is more strongly bound to the (HO)P 
unit at the expense of the CH 2 CHO— VO f bond, respectively 
(8 V2 p 2 : r((HO)P-OCHCH 2 ) = 1.691 A, r(O y V-OCHCH 2 ) = 
2.069 A; 8 V3P : r((HO)P-OCHCH 2 ) = 1.708 A, r(O f V- 
OCHCH 2 ) = 2.071 A). The relative energies of the respec- 
tive transition structures TS8-9 are however comparable 
(within a range of 12 kJmol" 1 ) when related to the entrance 
channel, respectively. As to the homonuclear phosphorus 
cluster [P 4 O 10 ]' + , intermediates 7 and 9 are linked via a 
single transition state; here, the second hydrogen atom is 
transferred directly to a terminal P=0 unit (see the Support- 
ing Information); however, this process is 67 kJmol" 1 higher 
in energy compared to the generation of the HAT product. 
For all vanadium-containing clusters, the spin distribution 
remains by and large unchanged along the reaction se- 
quence 7^TS7-8^8^TS8-9^9, carrying most of the spin 
density by a vanadium atom which is remote from the reac- 
tive sites. From 9, the weakly bound acetylene ligand can 
easily be liberated to yield the experimentally observed oxi- 
dative-dehydrogenation (ODH) product [V 2 P 2 O 10 H 2 ]' + 
(10). [26] As already mentioned above, an alternative reaction 
pathway has been identified for the ODH-process of the 
[V 2 P 2 O 10 ]' + cluster in which a structural isomer of 10 V2 p 2 , 
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Figure 3. Alternative mechanism for the ODH-reaction of [V 2 P 2 O 10 ]' + 
with C 2 H 4 , calculated at the B3LYP/aug-cc-pVTZ//B3LYP/TZVP level of 
theory (green V, yellow P, red O, gray C, white H). The electronic ener- 
gies and relative Gibbs free energies (in parenthesis) are given in 
kJmol 1 and corrected for unsealed zero-point energy contributions. The 
blue line shows the pathway from Figure 2. 



that is, 10y 2 p 2 -2 in Figure 3, is generated bearing two P— OH 
moieties instead of the P— OH/V— OH groups as in 10 V2 p 2 . 
The relative energy (free energy) of 10 V2 p 2 -2 is 153 
(150) kJmol" 1 lower compared to the one of 10 V2 p 2 ; however, 
the corresponding reaction mechanism generating the 
former product ion (Figure 3) is kinetically more energy de- 
manding than the HAT exit channel as well as the ODH 
pathway shown in Figure 2, that is, the hydrocarbon migra- 
tion from the terminal- to the bridging oxygen atom (TS7- 
8 V2 p-2) is 8 (42)kJmol" 1 and 36 (31) kJmol" 1 higher in 
(free) energy compared to the HAT product 6 V2 p 2 and to 
TS7-8y 2 p 2 -2, respectively. Alternatively, starting from 5 V2 p 2 
with a direct rebound of the C 2 fragment to the bridging 
oxygen atom of the P-O-P unit via 5 V2 p 2 ^8 V2 p 2 -2 might be 
possible but also here the PES is rather flat and the respec- 
tive transition structure could not be located. Intermediate 
8y 2 p 2 -2 corresponds to an open-cage cluster with a planar 
P0 3 unit generated by the cleavage of the internal P-O-P 
bond. Next, the transfer of a second hydrogen atom com- 
pletes the ODH reaction, and the process is accompanied 
with the regeneration of the closed-cage structure 9 V2 p 2 -2; 
the intrinsic barrier 8 V2 p 2 -2— >TS8-9 V2 p -2 is about 40 klmol -1 
lower than the related barrier 8 V2 p 2 -^TS8-9 V2 p 2 shown in 
Figure 2. 

As to the well-studied homonuclear clusters [V 4 O 10 ]' + and 
[P 4 O 10 ]' + , the reactivity patterns observed in the reaction 
with ethene are consistent with the relative energies of the 



competing pathways calculated in this study. For [V 4 O 10 ]' + , 
OAT is kinetically and thermochemically favored; the same 
holds true for the HAT pathway in the case of [P 4 O 10 ]' + for 
which no barrier exists along the reaction path. The latter 
process gains in importance in going from [V 4 O 10 ]' + via 
[V 3 PO 10 ]* + and [V 2 P 2 O 10 ]* + to [P 4 O 10 ]' + indicating that i) the 
hydrogen-atom affinity of the P— O' unit in [P 4 O 10 ]' + is 
higher compared to that of the V— O/ entity in [V 4 O 10 ]' + and 
that ii) the newly generated PO — H bond increases in 
strength with the number of phosphorus atoms being pres- 
ent in the cluster. In contrast, the reaction energies for an 
oxygen-atom transfer from [V 3 PO 10 ]' + , [V 2 P 2 O 10 ]' + , and 
[P 4 O 10 ]' + to C 2 H 4 are similar and OAT is most exothermic 
for [V 4 O 10 ]* + (Table 4). Note that the HAT and OAT reac- 
tion channels are different with respect to the redox chemis- 
try: While in the OAT process [V x P 4 _ x O 10 ] ,+ +C 2 H 4 ^ 
[V x P 4 _ x 0 9 ] ,+ +C 2 H 4 0 a reduction from +V to +IV takes 
place for V (x = 4) and P (x = 0, 2, 3), respectively, the 
formal oxidation states of the V and P atoms do not change 
in the course of HAT. 

The branching ratio of OAT versus HAT/ODH is mainly 
determined by the relative energies of the barriers TS2-3 
versus TS2-5 in Figure 2, respectively; both are lower in 
energy compared to the entrance channel for all clusters in- 
vestigated, Table 4. As mentioned above, the number of 
P atoms in the [V x P 4 _ x O 10 ] ,+ (x = 0, 2, 3) clusters has a negli- 
gible effect on the strength of the corresponding BDEs- 
(P— O f ) and thus affects the energies of the OAT products 
[V x P 4 _ x 0 9 ]" + (x = 0, 2, 3) only marginally. However, the rela- 
tive energies of the corresponding OAT transition structures 
TS2-3 increase discontinuously with the number of P atoms 
(Table 4). With respect to the HAT channel, the replace- 
ment of one vanadium atom by phosphor, that is, [V 4 O 10 ]' + 
versus [V 3 PO 10 ]' + , has only a minor effect on the relative en- 
ergies of the corresponding TS2-5 while the increasing ther- 
mochemical preference for HAT in going from [V 2 P 2 O 10 ]' + 
to [P 4 O 10 ]' + is even more pronounced when looking at the 
respective transition structures. As a result, the difference 
between the relative energies (free energies) of TS2-3 and 
TS2-5 amount to -65 (-58) klmol" 1 , -7 (0) kJmol 1 , and 0 
(7)kJmol" 1 , for the [V 4 O 10 ]' + /C 2 H 4 , [V 3 PO 10 ] ,+ /C 2 H 4 , and 
[V 2 P 2 O 10 ]' + /C 2 H 4 systems, respectively; for the [P 4 O 10 ]' + 
/C 2 H 4 -system there is no barrier for the HAT process. Thus, 
the calculations are consistent with the experimentally ob- 
served selectivities for OAT and HAT in the gas-phase reac- 
tions of [V 4 O 10 ]' + and [P 4 O 10 ]' + with C 2 H 4 , respectively. Al- 
though OAT in the reaction of [V 3 PO 10 ]* + and [V 2 P 2 O 10 ]' + 
with C 2 H 4 is exergonic and kinetically allowed at room tem- 
perature, the thermodynamic preference of the HAT and 
ODH products may account for the absence of OAT in the 
experiments. The non-observation of the second hydrogen 
transfer leading to ODH in the case of the homonuclear 
system [P 4 O 10 ]' + /C 2 H 4 , however, is in line with the calcula- 
tions since the corresponding transition structure TS8-9p 4 is 
much higher in energy compared to the HAT exit channel 
so that ODH cannot compete and is thus not observable 
under the experimental conditions. 
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Figure 4. Potential-energy surface for the reactions of [V 2 P 2 O 10 ] ,+ with C 2 H 6 , calculated at the B3LYP/aug-cc-pVTZ//B3LYP/TZVP level of theory 
(green V, yellow P, red O, gray C, white H). The electronic energies and relative Gibbs free energies (in parenthesis) are given in kJmol 1 and corrected 
for unsealed zero-point energy contributions. 



Ion/molecule reactions with C 2 H 6 : [V 2 P2O 10 ]' + has also been 
chosen as cluster prototype to demonstrate the reactions of 
the various cluster ions [V x P 4 _ x O 10 ]' + (x = 0, 2-4) with 
ethane, and the potential-energy profile is shown in Figure 4 
(Cartesian coordinates of all structures can be found in the 
Supporting Information and the relative electronic and free 
energies are summarized in Table 5). 

Similar to the reactions of these cluster-cation radicals 
with methane, coordination of C 2 H 6 to the clusters 
[V x P 4 _ x O 10 ]' + (x = 0, 2, 3) also results in a barrierless genera- 
tion of intermediate llp 4 , H V2 p 2 > anc * Hv 3 b respectively; once 



Table 5. Relative energies and free energies (parenthesis), given in kJmol , for the 
OAT, the HAT and the ODH channels calculated for the reactions of [V,P 4 x O 10 ]' + 
(x = 0, 2-4) with C 2 H 6 , at the B3LYP/aug-cc-pVTZ//B3LYP/TZVP level of theory. 



more, the P— O' unit constitutes the reactive site for the hy- 
drogen-atom transfer; in contrast to the reactions with 
ethene, no encounter complex has been located on the po- 
tential-energy surface. With respect to the homonuclear 
[V 4 O 10 ]' + cluster, the corresponding intermediate 11 V4 could 
not be located as a local minimum on the PES; instead, the 
hydrocarbon fragment migrates in the course of the HAT re- 
action around the terminal oxygen atom and rebinds to it 
thus forming directly intermediate 12 V4 bearing an intact 
ethanol ligand. Quite recently, the formation of CH 3 OH has 
also been described as a combination of HAT and CH 3 ' re- 
bound in the reaction of [Ai 2 0 3 ]' + with CH 4 under 
thermal conditions; in addition, this system also 
gives rise to the direct conversion of methane to 
formaldehyde. [27] 
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For the phosphorus-containing clusters, llp 4 , 
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and ll V3 p serve as common intermediates to 



[a] Directly linking intermediates 12 P4 and 16p 4 . 



branch out in three different reaction paths: first, 
the weakly-bound ethyl group can be eliminated 
yielding the cationic HAT product 6. As an alterna- 
tive, rebound of the ethyl group to the newly gener- 
ated OH group may occur to form C 2 H 5 OH which 
upon liberation of the alcohol gives rise to the OAT 
product 4. Here too, the transition state could not 
be located, the smooth PES in this region and the 
non-existing barrier for the 11 V4 ^12 V4 process indi- 
cate however a low activation energy for the re- 
bound step. The C— O bond lengths in 12y 4 , 12 V3 p> 
and 12 V2 p 2 increase with the number of phosphorus 
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atoms (1.517 A, 1.567 A, and 1.588 A, respectively) and all 
three are elongated compared to free ethanol (1.430 A); the 
C-O bond length in 12p 4 amounts to 1.559 A. While 12y 4 , 
12 V3 fs and 12 V2 p 2 correspond to closed-cage structures, a dis- 
torted open-cage structure was found for 12p 4 as the local 
minimum with one (C 2 H 5 OH)PO— P being broken, thus 
forming a bend P0 3 moiety with a threefold coordinated 
phosphorus atom at which 50 % of the spin is located. 

The third reaction path starting from phosphorus-contain- 
ing intermediates 12 comprises two variants of ODH reac- 
tions. In the first scenario, the open-cage structure 13 can be 
formed from intermediate 12 via TS12-13; this transition 
state corresponds to the second HAT from the terminal 
methyl group to the bridging P— O b unit thus leading to a 
(O fe ) 2 P(OH) 2 moiety to which the newly generated ethene 
molecule is only weakly coordinated. Regarding the elec- 
tronic structures of the vanadium-containing clusters, the 
spin density is located at one of the vanadium atoms 
throughout the whole reaction sequence. In contrast, the 
spin in TS12-13p 4 is transferred from the bent P0 3 unit to an 
adjacent phosphorus atom which is accompanied with a 
second P-O bond scission (r(P-O) = 1.607 A and 2.943 A in 
12p 4 and TS12-13p 4 , respectively), thus resulting in a planar 
P0 3 unit possessing two terminal oxygen atoms. The latter 
process requires 138kJmol" 1 activation energy and is thus 
23 kJmol -1 higher in energy compared to the entrance chan- 
nel. In all transition structures, one hydrogen atom of the 
terminal CH 3 group is transferred to the oxygen-atom and 
the respective C— H bond is elongated (r(C— H) 1.125 A for 
TS12-13p 4 , 1.118 A for TS12-13 V2 p 2 , 1.130 A for TS12-13 V3 p 
and 1.191 A for TS12-13 V4 , respectively). Liberation of C 2 H 4 
yields the product ion 14. 



In the second ODH-mechanism, the reaction sequence 
12^TS12 15^15 ^TS15 16^16^10 parallels the pathway 
for the oxidative dehydrogenation of ethene as discussed 
above (Figure 2). Briefly, for all vanadium-containing clus- 
ters the hydrocarbon fragment is first transferred to a bridg- 
ing oxygen atom via the isomerization sequence 12^TS12- 
15— >15; the barrier of this process is almost isoenergetic to 
the alternative transition state TS12-13 for [V 2 P 2 O 10 ]' + 
(+2kJmol" 1 ) and [V 3 PO 10 ]' + (+1 kJmol" 1 ), respectively, 
while it is distinctly lower for [V 4 O 10 ]' + (—28 kJmol" 1 ). Sub- 
sequently, the second HAT process occurs from intermedi- 
ate 15 via TS15-16 in which the hydrogen atom is transfer- 
red to the adjacent V— O t unit. Again, the homonuclear 
phosphorus cluster represents an exception in that the 
second C— H bond activation occurs directly from 12p 4 at 
one of the terminal P— O t units, without a prior migration of 
the hydrocarbon to a bridging oxygen atom. From 16, the 
ethene molecule is easily liberated, leading to the formation 
of the thermochemically preferred closed-cage product 10 
(Table 5). 

Similar to the [V 2 P 2 O 10 ]' + /C 2 H 4 couple, an alternative 
pathway for the oxidative dehydrogenation in the reaction 
of [V 2 P 2 O 10 ]' + with C 2 H 6 has been found in which the more 
stable product ion 10 V2 p 2 -2 (Figure 5) is generated bearing 
two P— OH units instead of one V— OH and one P— 
OH group. However, this process is kinetically more de- 
manding compared to the two ODH pathways described 
above. Starting from intermediate 12 V2 p 2 , the cage structure 
is opened via TS12-17 V2 p 2 , before the transfer of the ethyl 
radical to the newly formed P— O t unit via TS17-18 V2 p 2 takes 
place; the latter corresponds to the transition structure with 
the highest relative energy (— 182 kJmol" 1 ). The second 




Figure 5. Alternative mechanism for the ODH-reaction of [V 2 P 2 O 10 ]' + with C 2 H 6 , calculated at the B3LYP/aug-cc-pVTZ//B3LYP/TZVP level of theory 
(green V, yellow P, red O, gray C, white H). The electronic energies and relative Gibbs free energies (in parenthesis) are given in kJmol 1 and corrected 
for unsealed zero-point energy contributions. The blue lines show the initial pathways from Figure 4. 
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C— H bond activation of the oxidative dehydrogenation con- 
comitant with the regeneration of the closed-cage structure 
19y 2 p 2 occurs in the sequence 18v 2 p 2 -^TS18-19 V2 p 2 ^19 V2 p 2 . Fi- 
nally, elimination of the weakly-bound ethene brings about 
formation of 10 V2 p 2 -2; the overall process is associated with 
a significant thermochemical driving force of —349 kJmol" 1 . 

While the reaction channels discussed represent the ener- 
getically favored pathways for ODH, many other isomeric 
intermediates are possible as well, depending on the coordi- 
nation of the hydrocarbon, the hydrogen atom, etc. Howev- 
er, the rate-determining activation barriers of the two ODH 
pathways determined for the [V 3 PO 10 ]' + /C 2 H 6 couple are 
TS12-13 V3 p and TS12-15 V3 p. While these transition structures 
with relative energies of —169 and — 168kJmol" 1 , respec- 
tively, are almost isoenergetic, the associated product ions 
are quite different energetically with the formation of 10 V3 p 
being significantly more exothermic, that is, 138 kJmol 1 , 
compared to 14 V3 p. For the [V 4 O 10 ]' + /C2H 6 system, 10 V4 is 
not only thermochemically favored by approximately 
154kJmol" 1 , but also the respective transition state TS12- 
15y 4 is 28 kJmol" 1 lower in energy than the alternative tran- 
sition state TS12-13y 4 leading to 14 V4 . Moreover, formation 
of 10y 4 has been demonstrated in the reaction of [V 4 O 10 ]' + 
with C 3 H 8 in an elegant combined experimental/computa- 
tional study including infrared spectroscopic characteriza- 
tion of the ionic species. [18d] 

In line with our experiments (Table 1), formation of the 
OAT product 4 has been calculated to be endothermic by 
14, 12, and 16kJmol" 1 for heteronuclear [V 3 PO 10 ]' + and 
[V 2 P 2 O 10 ]' + and homonuclear [P 4 O 10 ]' + , respectively. In con- 
trast, OAT is exothermic by -103 kJmol" 1 and thus accessi- 
ble for the homonuclear vanadium [V4O 10 ]' + /C 2 H 6 system 
under thermal conditions. However, the ODH reaction 
channel for this system is even more exothermic by 
183 kJmol" 1 , and the highest barrier of this pathway, that is, 
TS12-15y 4 , has been calculated to be about 83 kJmol" 1 lower 
in energy compared to the exit channel of the oxygen-atom 
transfer. Accordingly, consistent with theory and our experi- 
mental results, [19] oxidative dehydrogenation corresponds to 
the main reaction channel in the reaction of [V 4 O 10 ]' + with 
C 2 H 6 . 

Conclusion 



more reactive with respect to homolytic C— H bond cleavage 
to form the closed-shell HAT product [P 4 0 9 (OH)] + . Fur- 
ther, while the first C— H bond scission is initiated by an 
active P— O' site, the presence of a vanadium atom in the 
cluster is essential to enable the second hydrogen-atom 
transfer from the hydrocarbon, that is, oxidative dehydro- 
genation. Obviously, while the presence of vanadium is cru- 
cial a direct interaction of the redox-active transition metal 
with the hydrocarbon is not required to bring about ODH 
reactivity under thermal conditions. Thus, cooperative ef- 
fects of V and P atoms in the cluster are operative in the oc- 
currence of ODH; neither [V 4 O 10 ]' + nor [P 4 O 10 ]' + alone pos- 
sess the electronic properties necessary to bring about the 
dehydrogenation of ethene. 

The redox activity of the P and V atoms, that is, accessibil- 
ity of the oxidation state +IV, is also of importance in the 
reaction of the ionic cluster with ethane. Due to the prefer- 
red oxidation state +V, [P 4 O 10 ]' + is much more reactive in 
terms of homolytic C— H bond activation during which the 
formal oxidation state does not change, in contrast to OAT 
in which a phosphorus atom is necessarily reduced. On the 
other hand, vanadium can be reduced from + V to +IV and 
can therefore easier stabilize the radical site formed in the 
course of reductive processes. Furthermore, the +V— >+IV 
reduction of vanadium is associated with rather small struc- 
tural rearrangements, causing only a insignificant elongation 
of the internal V— O bonds of the cluster; accordingly, 
[V 4 O 10 ]' + favors the formation of the open-shell products 
[V 4 0 9 ]' + and [V 4 O 10 H 2 ]' + . Recalling the initial questions, 
some crucial aspects about the intrinsic features of the 
mixed-systems with respect to C— H bond activation have 
been revealed based on these characteristic features of phos- 
phorus and vanadium. Thus, the combination of both ele- 
ments gives rise to new product distributions and illustrates 
the cooperative effects between different metals and non- 
metals in complex oxo-frameworks. Although the charge of 
the gas-phase clusters can affect the energetics of the inves- 
tigated processes compared to the neutral systems, [28] the in- 
vestigation of model systems can help to elucidate reaction 
mechanisms and give hints with respect to the structure of 
reactive intermediates in heterogeneous catalysis ; [4 29] these 
insights should be considered relevant for real catalytic sys- 
tems. 



We have analyzed the reactions of various hetero- and ho- 
monuclear phosphorus-vanadium oxygen-cluster ions 
[V x P 4 _ x O 10 ] ,+ (x = 0, 2-4) with ethene and ethane using 
DFT. These calculations, conducted at the B3LYP/aug-cc- 
pVTZ//B3LYP/TZVP level of theory, reveal some interest- 
ing aspects concerning the underlying reaction mechanisms 
which are operative and which explain the puzzling reactivi- 
ty patterns observed in the gas phase under ambient condi- 
tions. While [V 4 O 10 ]* + reacts with ethene exclusively in 
terms of oxygen-atom transfer to the hydrocarbon, the pres- 
ence of already one phosphorus atom inhibits this pathway 
due to the strong P— O* bond. In contrast, [P 4 O 10 ]' + is much 



Computational Methods 

All calculations were performed by using the Gaussian09 package. [30] Ge- 
ometries were optimized at the unrestricted UB3LYP level of theory [31] 
with the triple-^ plus polarization basis sets TZVP. [32] B3LYP was shown 
previously to describe vanadium-oxide clusters in good agreement with 
available experimental data as well as quantum chemical methods that 
explicitly include electron correlation. [33] For example, B3LYP repro- 
duced experimental and CCSD(T) results for the dissociation of [V0 2 ] + 
into [VO] + and 1 / 2 0 2 ; [29cJ good agreement between multi-reference calcu- 
lations and B3LYP results has also been found for the molecular struc- 
tures of V 2 0 4 and the relative energies of its open shell singlet and trip- 
let. [34] Thus, this method proved reliable in numerous studies of cationic 
oxo-cluster ions, containing vanadium and/or phosphorus. [9d l8c ' 20ab ' 29c ' 33 " 35] 
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Vibrational frequency analyses have been carried out at the same level of 
theory to characterize the nature of stationary points as minima or transi- 
tion structures, and to derive the zero-point energy corrections (ZPE). 
Further, the energy of all stationary points was computed also by single- 
point calculations, using the larger, correlation-consistent basis set aug- 
cc-pVTZ of Dunning and coworkers, to derive more accurate values for 
the apparent activation barriers. [36] All relative energies (corrected for 
ZPE contributions) and Gibbs free energies (at STP) are reported in 
kJmol -1 . Intrinsic reaction coordinate (IRC) [37] calculations or manual 
displacement along the reaction trajectory of the imaginary frequency 
were performed to link transition-state structures with the respective in- 
termediates. 
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